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Abstract 

We present a theory for the dynamical evolution of a quantum system coupled 
to a complex many-body intrinsic system/environment. By modelling the 
intrinsic many-body system with parametric random matrices, we study the 
types of effective stochastic models which emerge from random matrix theory. 
Using the Feynman- Vernon path integral formalism, we derive the influence 
functional and obtain either analytical or numerical solutions for the time 
evolution of the entire quantum system. We discuss thoroughly the structure 
of the solutions for some representative cases and make connections to well 
known limiting results, particularly to Brownian motion, Kramers classical 
limit and the Caldeira-Leggett approach. 



PACS numbers: 05.60+w, 05.40.+j, 24.10.Cn, 24.60.-k 
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I. INTRODUCTION 



Quantum dissipation is a problem with such a long history and such a multitude of results 
that even a cursory review will not make justice to the numerous contributions of a large 
number of authors over several decades. We recommend the reader to consult Refs. flj-pUTT 



and the references therein. In spite of this impressive effort, many people do not consider the 
problem of quantum dissipation a solved issue, and its character and microscopic origin still 
call for the attention of a large community across many, if not all, subfields of physics. The 
present paper represents a continuation of our effort to understand the character of energy 
flow between the slow degrees of freedom and the intrinsic degrees of freedom in many- 
body systems. Our initial motivation was to understand the "irreversible" time evolution of 
the large amplitude nuclear collective motion p|-|TT)| and for that reason we have adopted a 
traditionally nuclear physics approach. The large body of evidence, both experimental and 
theoretical, suggests that many Fermion systems can be described reasonably well within the 
framework of a random matrix formalism. Another way of saying the same thing is that a 
many Fermion system is predominantly a quantum chaotic system and thus a random matrix 
approach is a natural approach. At the same time, both theory and experiment strongly 
suggest that there are some degrees of freedom, which are not chaotic and usually are referred 
to as collective or shape. However these collective degrees of freedom are coupled with a 
large number of non-collective degrees of freedom and as a result a rather generic situation 
results: a relatively small quantum system in contact with an "environment" . Even though 
the whole system is finite and in a strict sense there is no irreversible behaviour in this case, 
for all practical purposes the time evolution of the collective or slow quantum system has 
the character of quantum dissipative dynamics. In our formulation, we will neglect other 
physical mechanisms which lead to dissipative contributions, such as particle evaporation 
or coupling to electromagnetic fields. This can be viewed as limiting the present study to 
shorter time scales. Such phenomena can in principle be introduced into the formalism, but 
we shall not attempt it here. We shall not either try here to defend the legitimacy of such a 
terminology (dissipation), the issue at stake however is unquestionably sound. The reader 
will recognize easily that the problem we address here is typical and under different guises 
appears in many subfields of physics. 

As an introduction, let us consider for the moment that a certain simple system interacts 
with some relatively large (but finite) many-body system. The question is: can one describe 
the dynamical behaviour of the simple system using for example an equation of the form 

„ d 2 X dU(X) t , dX £lv , „ 

as in the case of a Brownian particle, if in the absence of the interaction the Hamiltonian of 
the system is 

H (X,P) = ^ + U(X) (1.2) 

and where 7 is a friction coefficient and f(X,t) is a Langevin-like force? The force f(X,t) 
can in principle depend not only on time but on position as well, and in this way one can 
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describe a large variety of physical situations, ranging from diffusion to localization ||11|| . If 
one were to start from a description of the entire system (reservoir plus simple system) with 
a Hamiltonian 

H(X, P, x,p) = H (X, P) + H X (X, x,p), (1.3) 

where Hi(X,x,p) describes the reservoir and its interaction with our system, under what 
circumstances could one derive an equation of motion like Eq. ( [Lip ? Moreover, does the 
fluctuating force have Gaussian character or not? Since our approach is a fully quantum 
mechanical one, we shall be able to answer another important question, on the specific 
role played by the quantum effects. As we shall show, the classical results are recovered 
when the appropriate limits are taken, however, the classical limit is not at all a trivial 
one nor is it reached in a simple fashion. If one were to take the approach of postulating 
that an equation of the type Eq. (|1.1|) governs the dynamics and assume for example 
that the fluctuation properties of f(X, t) are Gaussian in character, then the entire powerful 



apparatus developed for Brownian motion could be then invoked [[LT| . But such an approach 
would leave unanswered the main questions of whether one can describe in this way collective 
degrees of freedom in a finite closed system, like atomic nuclei. 

We address this problem using a well known approach based on the double path integral 
formulation of Feynman and Vernon []]]. Our original input is in the functional form of the 
influence functional, which arises from a parametric random matrix description of the "envi- 
ronment". This has been attempted earlier in Ref. ||. The functional form for the influence 
functional we determine is qualitatively different from the popular Caldeira-Leggett type 
derived by Feynman and Vernon The parameters which define the influence functional 
have a rather simple and intuitive meaning from a microscopic point of view and we refer 
the reader to earlier publications for details and discussions []5| ITOfl • It comes as no surprise 
than that under such circumstances the dynamical evolution of a quantum dissipative sys- 
tem in our case has new features as well, as we shall amply exemplify in the body of the 
paper. Here we restrict our attention to the Markovian limit only and we hope to address 
the important problem of memory effects in the future. In spite of its physical restrictions 
(high temperature limit for the intrinsic system) this limit shows already the qualitative 
differences with the previously known approaches. 

The paper is organized as follows. In Section II we discuss the time evolution equation for 
the density matrix of the "slow" quantum system coupled to a complex many-body system. 
In Section III we show that, at high temperatures, the evolution equation for the density 
matrix can be brought to the Kramers form, when the classical limit is taken. Sections IV- 
VII discuss exact solutions to the evolution equations for certain potentials. Section VIII 
is devoted to the study of the tunneling problem. A short summary and discussion of the 
results is given in the final Section. 



II. EVOLUTION EQUATION FOR THE DENSITY MATRIX 

In this section we discuss the description of the internal degrees of freedom (or complex 
environment) through parametric random matrix theory, and derive the equation of motion 
for the density matrix of the slow degrees of freedom by integrating over the internal states. 
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A. Random Matrix Model 



The basic assumption concerning the intrinsic states is that there are no governing con- 
stants of the motion, so that the dynamics is chaotic. This has been seen to be the general 
situation in studies of many-body systems, from nuclei to molecules, so it is reasonable to 
approach the modelling of these degrees of freedom with random matrices, suitably tailored 
to the problem. The generic form of the Hamiltonian governing the dynamics of a "slow" 
quantum system coupled to a complex environment is described as follows 

H(X, x) = H (X) + H^X, x) = ^ + U(X) + H X (X, x). (2.1) 

We shall often refer to X as "shape" variables, since in large amplitude collective nuclear 
motion it represents the collective coordinates which characterize the nuclear mean field. 

The part of the total Hamiltonian Eq. ( [2.1Q which depends on the intrinsic coordinates 
H\{X, x) is defined as a matrix, whose matrix elements depend parametrically on the "slow" 
coordinate X 

[H 1 (X)) ij = [h ) ii +[h 1 {X)] ij . (2.2) 

ho is taken to be diagonal and defines the average density of states, with (k\h \l) = [ho]ki = 
EkSki- We refer in the main text to these eigenstates as "typical states" of the intrinsic system 
with an energy e. One can think of ho as a Hamiltonian describing a "bath" or a "reservoir" 
and of hi(X) as a Hamiltonian describing the interaction between the "bath/reservoir" and 
the "slow" system. Whereas in statistical physics the interaction between the thermostat 
and the system under consideration is assumed to be negligible, we shall not make such an 
approximation here. As a matter of fact, for the physically interesting situations we envision, 
this coupling term can be large. This fact alone leads to significant differences of various 
distributions when compared with the corresponding results of traditional approaches. 

For an intrinsic subsystem with a large number of degrees of freedom, the average density 
of states, 

p(e) = Tr5(H!(X)-e), (2.3) 

for each given shape X increases sharply with energy. The overline denotes here a procedure 
for extracting the smooth part of p(e) as a function of energy, which amounts essentially 
to an ensemble average, to be introduced below. For a many Fermion system, p(e) has 
a roughly exponential behaviour. Recall that \np(e) is approximately proportional to the 
thermodynamic entropy of the intrinsic system, which is an extensive quantity. The fact that 
the average density of states for the intrinsic subsystem has such a behavior is a key element 
of the entire approach. This is equivalent to stating that the intrinsic subsystem has a large 
heat capacity and thus can play the role of a "reservoir" , although not necessarily ideal. In 
principle p(e) can be A-dependent as well, but we shall ignore this aspect here. Without 
an A-dependence of the average density of states, mechanical work cannot be performed on 
or by the model environment we study here, and only heat exchange is allowed. 

In Refs. [|6|-|T0|j we discussed the reasons why one chooses this specific form of the Hamil- 
tonian. In the basis of the eigenstates of h Q , we define hi(X) as a parameter dependent, 
N x N real Gaussian random matrix, which is completely specified by its first two moments 
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MX)] 



kl 



M*)kM^)W = [SikSn + SuS^GijiX - Y). 



(2.4) 



The overline stands for statistical averages over the ensemble of random Gaussian matrices 
from the Gaussian Orthogonal Ensemble (GOE) fT2"|l . Qij(X — Y) can be taken as a "bell 
shaped" correlation function with a characteristic width X Q , or, in some physically interest- 
ing cases even periodic, with period oc Xq. The fact that this correlator is "translational 
invariant" is not a crucial limitation, and a general form can be adopted without any sig- 
nificant changes in the formalism. We limit our analysis to the GOE case only for the sake 
of simplicity of the argument, as any other Gaussian ensemble can be treated in a similar 
manner. The dependence on i,j allows for the description of banded matrices, where an 
effective number of states Nq < N are coupled by hi(X). It is convenient to use an explicit 
parameterization, which incorporates the average density of states and the bandwidth of the 
statistical fluctuations explicitly |5|]: 



QiAx) 



2^^p(e i )p(e 3 



exp 



2 Kg 



(2.5) 



Here G(x) = G(—x) = G*(x) < 1, G(0) = 1, and the spreading width I^, k (linked with 
the effective band width iVo « K, p(e)) and X Q are characteristic of the intrinsic system. 
Even though it is not necessary, in this paper we shall consider a particular from for G(x), 
namely G(x) = exp(— x 2 /2). 

To complete the tailoring of the random matrices, we require a realistic average density 
of states for the reservoir. It is reasonable to assume that in a suitable energy interval, p 
has the behaviour [Hi 



p(e) = p exp (/fe), 



(2.6) 



(3 can thus be interpreted as the thermodynamic temperature of the intrinsic system. As 
(3 is independent of the internal excitation energy e, this particular type of intrinsic quan- 
tum system plays the role of a "perfect" thermal quantum reservoir. We note in passing 
that even though the temperature of the "reservoir" remains constant throughout the entire 
dynamical evolution of the whole system, one should not conclude from this that the "reser- 
voir" is in thermal equilibrium. As we have shown explicitly in Refs. |||| the population of 
various energy levels of a uniformly driven "reservoir" is far from an equilibrium Boltzmann 
distribution. The parameters of the present construction are summarized in Table I. 



B. Influence Functional 



The quantum description of our coupled system will be treated through the path integral 
construction of the density matrix. According to Feynman and Vernon |TJ, one can write 
the following double path integral representation for the density matrix of the entire system 
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K(X,x,Y,y,t)= dX dY iP(X )r(Y ) 



X{t)=X 
X(0)=X 



VX(t) 



Y(t)=Y 
Y(0)=Y 



VY(t) 



xexp\-[S (X(t))-S (Y(t))} 



x (x|Texp 



4 f dt'H^Xit')) 
n Jo 



T a exp 



ifjfH^Yit")) 



(2.7) 



where T and T a represent the time ordering and time anti-ordering operators respectively. 
In this representation, we have used a particular form for the initial state wave function, 



*(X,x) = ip(X)</>(x). 



(2.8) 



Other choices are equally possible, such as an initial density matrix. By introducing the 
influence functional £(X(t),Y(t),t) 



C(X(t),Y(t),t) 



T a exp 



dfH x {Y{t")) 



jjTexp 



h Jo 



dt'H^Xit')) 



(2.9) 



one readily obtains the following double path integral representation for the density matrix 
for the "slow" subsystem 



p(X,Y,t)= dX dY^{X )i,*{Yo) 



X(t)=X 
X(0)=X 



VX(t) 



Y(t)=Y 
Y{0)=Y 



VY(t) 



xexp\- [S (X(t)) - S (Y(t))} \ C(X(t), Y(t),t). 



(2.10) 



The formulation of the problem through a path integral representation serves only as a very 
convenient vehicle to obtain an evolution equation for the density matrix p(X,Y,t). 



C. Evolution Equations 



The evolution equation for the influence functional C(X(t),Y(t),t) has been solved in 
Ref. [j§|-|lOH for the case N — ► oo and the case when the "temperature" of the reservoir is 
infinite. In Appendices A-D we compute the first order correction in (3 to the influence 
functional in the adiabatic limit, when the characteristic time scale of the "reservoir" K/kq 
is significantly shorter than the characteristic time scale of the slow system for which we 
derive the dynamical evolution equation. We thus obtain for the influence functional 



C(X(t),Y(t)) = exp{^- f da 
{ n Jo 



G 



'X(s) - Y(s 
. An 



x exp 



4A Jo 



ds[X(s) +Y(s)]G' 




(2.11) 



where G'(x) = dG(x)/dx. The physical significance of all other quantities entering this ex- 
pression has been explained and discussed at length in Refs. 0-0, and is briefly summarized 
in Table I. It is worth noting that the functional form of the influence functional derived by us 
is different from the Caldeira-Leggett form H, which is a quadratic expression in X(t) and 
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Y(t). If we were to use only the first term in a Taylor expansion of G[(X(s) — Y{s))/Xq\ — 1, 
we would obtain an expression similar to Caldeira-Leggett form for the influence functional. 
However, the present form of the influence functional leads in the classical limit to a velocity 
dependent frictional force, see Ref. 

By combining the double path integral representation for the density matrix p(X, Y, t) 
with the above expression for the influence functional in the adiabatic approximation one 
readily obtains that the density matrix satisfies the following Schrodinger like equation (for 
similar examples see Refs. MM) 



G 



(X-Y 

v x 




4X M V X 

with an arbitrary initial condition 

p(X,Y,0) = po(X,Y). (2.13) 

This equation is the central object of our study and the remaining of the paper is devoted to 
determining various limiting regimes and the character of its solutions upon varying j3, 
Xo and H. (In this paper we have already taken the limit kq — > oo.) At first glance the reader 
might get the impression that the slow subsystem we consider here is characterized by one 
degree of freedom only. As a simple analysis will show however that the slow subsystem can 
have an arbitrary number of degrees of freedom and most of the formulas we shall present 
are equally valid in this case. 

This evolution equation is somewhat peculiar in certain aspects. It is obvious that in 
the absence of the coupling to the "reservoir", it describes a purely quantum evolution of 



the "slow" subsystem. Eq. (|2.12|) has been derived from a purely quantum description of 



the entire system, by performing the expansion in (3. In Eq. Q2.12 ) however, the inverse 



temperature enters only in the combination r t h = /3fi, which can be interpreted as a thermal 
time (analogous to the thermal de Broglie wave length). Thus the expansion in the inverse 
temperature f3 is at the same time an expansion in H. Our limitation to the zeroth and first 
order terms in (5U for the coupling between the "slow" subsystem and the "reservoir" can 
consequently be interpreted as a semiclassical approximation. 

It has been argued by Diosi || that for the case of a Caldeira-Leggett correlator (G(x) = 
1 — x 2 /2) the similar high temperature limit of the evolution equation requires the retention 
of the next order term in (3 in order to bring the corresponding approximate evolution 
equation to a Lindblad form fT3"| , which maintains positivity of the density distribution 
for any physically acceptable initial conditions. If these higher order terms in (3 are not 
introduced, such an equation (as Eq. ( j2.12| )) cannot be applied to an initial state which 



is narrower than the thermal de Broglie wave length, = j3 / M . This restriction 

to wave packets which are wider than the thermal de Broglie wave length is manifest in a 
somewhat different way as well. Let us compute the rate of change of the "total energy" of 
the "slow" subsystem, which can be defined naturally as follows: 

E (t) = Tr[H oP (t)]. (2.14) 
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Using the evolution equation Eq. (|2.12|) for the density matrix p(X, Y, t), after some straight- 
forward manipulations one obtains the following expression for the rate of change of the 
"total energy": 



dE (t) 
dt 



27 



T 
~2 



<P*> 
2M 



(2.15) 



where 7 is the friction coefficient in the small velocity limit to be introduced below, see 
Eq. (^). This rate has an apparently pure classical content. This is of course deceiving, 
as quantum effects are clearly retained in both Eqs. ( |2.12|) and ( |2.15| ), even though not 
entirely in Eq. fl2.15| ). However, since one has to assume that T > h 2 /ML 2 , where L is 
the characteristic spatial extension of the state, it is possible in Eq. (2.15) to replace the 
quantity (P 2 ) = (P) 2 + ((P 2 )) with simply (P) 2 . This renders Eq. (2.15) purely classical 
in character. One incurs a certain loss of accuracy and a small degree of inconsistency by 
proceeding in this manner, so it is better to leave Eq. (|2.15|) as is. 

Even though one can go beyond the first order in (3 and one can derive a more accurate 
evolution equation for the density matrix p(X,Y,t) in the high temperature limit, we shall 
not do that in this work, for the sake of simplicity of the presentation. We do not expect 
that such corrections will lead to a qualitatively new behaviour. 



D. Coordinate and Momentum Distributions and the Cumulants Expansion 

We will derive the time dependent solutions p(X, Y, t) in the following sections, and from 
that it will be useful to extract information concerning the behavior of coordinates and 
momenta. The most convenient way to do so is through the cumulant expansion. To define 
this, we start by introducing the new variables: 

r = X ~^ Y , s = X-Y. (2.16) 

Coordinate and momentum information can readily be extracted from the following rather 
atypical Fourier transform of the density matrix: 

f dk ( ikr\ ,. , . , 
p{r,s,t)= J — expf— (2-17) 



For either s = or k = 0, d(s,k,t) is the characteristic function [13] for the spatial or 
momentum distribution of a given quantum state respectively. For example, if we are in- 
terested in the spatial diffusion X, as measured by ((X 2 )), then we get from Eq. ( |2.17| ) and 
integration by parts: 

(X 2 ) = J dX'p(X', X', t)X' 2 = ~h 2 -^d(0, k, t) \ k=0 = ((X 2 )) + (X) 2 . (2.18) 
Similarly, in order to compute the average collective energy one needs (P 2 )'- 
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(P 2 ) = J dXdX'p(X,X',t)(X' \P 2 \X) (2.19) 



— / dXdX'dPp{X, X', t) exp 
27m J 



iP(X' - X) 



h 

= J dsdPexp d(s,0,t)P 2 = -h 2 ^d(s,0,t) \ s=0 = ((P 2 )) + (P> 2 . 

The quantities denoted ((••■)) are the cumulants of the distribution. 

From the definition of the function d(s, k, t) one defines the general coordinate and mo- 
mentum cumulant expansion as: 

lnd( S) M)| s=0 = $:-(j-j «0> (2.20) 

M«,M)k=o = £z7(i) M (2.21) 

n=l n - vm/ 

where ((r n )) and are the (time dependent) cumulants of the spatial and momentum 

distribution respectively. One can show that the zeroth order terms in both cumulant 
expansions vanish, which is consistent with the fact that the probability is conserved within 
the present formalism. 

A Gaussian process has only non-vanishing first and second cumulants. In general it is 
known from Marcienkiewics' [|17] or Pawula's |14| theorem, that for a probability distribution, 
one either has a Gaussian process with only the first two cumulants non-vanishing, or 
all cumulants are present. Furthermore, while there are some inequalities which relate 
cumulants of varying order, in most cases there is no restriction on their sign, which can be 
positive or negative. 



P 2 



E. Intrinsic States 



By "integrating" over the intrinsic variables one of course looses a lot of information and 
only some average behavior of the slow subsystem can be thus inferred. In principle one 
can recover some of this information in the following manner. As in Refs. one can 

introduce instead the generalized occupation number probabilities 



Af(X(t),Y(t),t,e) = (</>\T a exp 
x (e:|Texp 



dfH^Ytf')) 
' dt'H l (X(t')) 



(2.22) 



where \e) is a "typical" state of the intrinsic subsystem with energy e. One can derive 
now evolution equations for these quantities. With the help of these generalized occupation 
number probabilities one can introduce what was called in Refs. [§,|9| the characteristic 
functional 



N{X{t),Y{t),t,T) = J dep(s)Af(X(t),Y(t),t,e) exp 



(%ET 



(2.23) 
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where p(e) is the average level density for the intrinsic subsystem. From this expression, 
by performing an inverse Fourier transform one can reconstruct the generalized occupation 
probabilities Af(X(t),Y(t),t, e). In such a case there is a double path integral representation 
for the density matrix for the slow subsystem with the condition that the intrinsic subsystem 
has a well defined energy e, namely 



p(X,Y,t,e) = I dX dY 4j(X )r(Yo) [ VX{t) f VY(t) 

J Jx(0)=X JY{0)=Y o 

x exp{^[S (X(t)) - S (Y(t))]} p(e)Af(X(t),Y(t),t,e). (2.24) 



In principle this representation can be used to extract more detailed information concerning 
the time evolution of the coupled slow and complex intrinsic subsystems. We shall not 
attempt to do that however in this paper. 

III. CLASSICAL LIMIT: KRAMERS EQUATION 



It is interesting to explore the classical transport equation which emerges from Eq. ( |2 . 1 2 ) . 



The standard approach is to construct the Wigner transform f(Q } P, t) of the density matrix 
p(X, Y, t) as 

It is well known that while f(Q, P, t) can be interpreted as a classical probability distribution 
in a phase space (Q, P), it is only a quasi-probability since its sign can be positive or negative 
at a given phase space point, while its integral over any unit phase space cell (of size K) is 
positive semi-definite. We will further introduce the friction coefficient: 



^2 
L 



-' = fuxr (3 ' 2) 



As we will see, this definition will emerge naturally from our analysis of the dynamical 
evolution of the quantum systems. However, in taking the Wigner transform of p(X, Y, t), 
we must also take the classical limits of the quantities and X , which have an intrinsic 
quantum interpretation. While it is not entirely clear to us how to define such limits, we 
will see that the combination which appears in 7 has a natural classical interpretation. 

From the definition, one can readily compute the classical evolution equation, which is 
simply the Wigner transform of the right hand side of (|2.12 ): 



d t f{Q,P,t) = ^-J dRexp (-^-) d tP (Q + -R,Q- (3.3) 

- i^hx.G' (jp) d R } p (q + 1 -r,q- ±r, tj 



10 



For the contribution due to the potential energy U(X), one can use the Kramers-Moyal 
expansion 



dR ( iPR\ 

2^ exp r^rJ 



i (q + Ir\-u(q-~r 



(Q ■ -i/,'./; (3.D 



= U(Q) 2 - sin mP ,t) 

= d Q U(Q)d P f(Q,P,t) + o(K) 

where Oq and dp in the sine term act only on U and /, respectively. In the last line only 
terms to o(h) were retained. For the terms which depend upon the correlation function 
G(x), we consider a general expansion 

G(s)«l-y + ... (3.5) 

The terms beyond the quadratic ones in the expansion will come with higher powers of h 
into the evolution equation and are hence omitted. Integrating by parts and using the fact 
that p(X, Y, t) vanishes when X, Y — > ±oo, we have: 

df(Q,V,t) df(Q,V,t) 1 dU(Q)df(Q,V,t) 

dt dQ M OQ dV 1 ' ' 

\ d[Vf{Q^t)] T d 2 f(Q } V,t) } 
" 7 \ 8V M dV 2 J ' 

where T = 1/(3 is the temperature, and the velocity is V — Q — P/M. We have thus 
obtained Kramers equation flUR . 

We will derive transport coefficients below, and it is worth noting that taking classical 
limits is not straightforward. In quantum Brownian motion, we will extract a diffusion 
constant which is related to the friction 7 through the classical Einstein relation: 

_ 2X1 _ (37) 

We note here that our transport theory has a consistent classical limit for all of these 
transport coefficients only when they remain finite as Ti — > 0. This requires in turn that 
the parameters of our theory cannot remain constant as H —* 0, if we are to recover a well 
defined classical transport. 



IV. TIME DEPENDENT SOLUTIONS 

In the remaining part of this article we will consider solutions to the evolution equations 
(|2.12|) . For certain forms of the potential, one can readily obtain explicit solutions to the time 



evolution of the density matrix, while for others we solve the evolution equation numerically. 
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A. Exactly Solvable Limits 



For the class of potentials: 

f —FX Section V 

U(X) = < \Muj 2 X 2 Section VI . (4.1) 

[ -\M£l 2 X 2 Section VII 

the evolution equation for the density matrix can be written as a first order partial differential 



equation in two variables, which is readily solved by the method of characteristics [y| . 



B. Numerical solutions 

For other cases where one does not have analytic results, the method of solution can be 
chosen in many ways. For the double well potential, we have implemented two independent 
methods which agree. The first is a brute force technique in which we express the density 
matrix in the form 



p(X, Y,t) = J2 a nm (t)MX)(j>* n (Y) exp 



i{E m - E n )t 

h 



(4.2) 



where the 4>(X) are the eigenstates in the absence of coupling to the intrinsic degrees of 
freedom: 

H (X)<j> n (X) = E n <j> n (X). (4.3) 

The evolution equations for a mn (t) are solved numerically in a truncated dimensional 
basis, with m, n = 1, N max . Typically we will use N max = 31. 

Alternatively, we have solved Eq. (|2.12|) on a two dimensional grid (A% YJ). It is found 



that the time iteration is most stable using the alternating-direction implicit method JIB 



with five points interpolation for the space derivatives. As a check on the validity of the 
solution, the total probability is verified to be preserved at unity: 



dXp(X,X,t) = 1. (4.4) 
The latter method allows one to integrate much farther in time. 

V. LINEAR POTENTIAL AND QUANTUM BROWNIAN MOTION 

The classical picture of a Brownian particle in a constant force field F and interacting 
with a heat bath is described by the Langevin equation for the velocity: 

v + lv -L = f{ t) (5.1) 

Here 7 is the friction coefficient, and f(t) is Gaussian white noise. In the long time limit the 
particle energy equilibrates with ((p 2 )) = TM and there is a finite drift velocity = F/jM. 
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In this section we consider the dynamics of a quantum particle in a constant force field 
interacting with our random matrix bath and contrast it to this classical limit. The results 
discussed here should be contrasted with those for the case /3 = 0, where we have found that 
the quantum dynamics is similar to turbulent diffusion PIJTO . 



A. Exact solution 

Let us consider the case when there is a linear potential acting on the slow variables 



H (X) 



2M 



:& 



x 



FX. 



(5.2) 



One can consider the case of a time dependent linear potential as well, i.e. U (X, t) = 
—F(t)X, with only very slight modifications of the ensuing formulas. The equation for the 
density matrix now becomes 



inut + -T70 r O s - 

M 2X M 



O' ( ^ ) 0, 



p(r, s, t) = {-Fs + iT l [G(— ) - l]}p(r, s, t) (5.3) 



with the initial condition p(r, s, 0) = po( r , s)- It is not necessary to consider a pure state 
as an initial state, and we allow for any general initial density matrix. The mixed partial 
derivative can be removed by passing to the Fourier transformed equation (see Section II. D) 
for d(s, k,t). This satisfies the equation 



d t + 



A _ ElLa (A. 

M 2MX \X 



d s } d(s,k,t) 



iFs r l 



G 



X 



d(s,k,t). (5.4) 



Using the method of characteristics for wave equations |15[] one can find the solution in 
parametric form: 



d(s,k,t) = d (S(t), k) exp < J dt' 



iFS(t - 1') 

h + T 



G 



'Sjt-t'Y 

Xa 



- 1 



p(r,s,t) 



dk 
2nh 



d (S(t),k) 



x exp ^ I dt' 



iFS(t - o 



G 



'S(t - 1') 

x n 



(5.5) 



(5.6) 



where the time dependent function S(t — if) is the solution of the auxiliary equation 

dS{r) 



M 2MXn 



Xn 



(5.7) 



Here d(s, k, t)\ t=0 = do(s, k) is the initial distribution, which is just the Fourier transform of 
Po(r,s). The coordinate s appears in this solution as an initial condition on S: S(0) = s. 
These equations define the flow of the density matrix in time. For an arbitrary if the function 
S(t — if) also satisfies the homogeneous equation 
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d t + 



k pT^h ( s 



M 2X M 



Xr 



d s \ S(t - t') = 0. 



(5.8) 



Again, S(t — t') depends on s through the initial condition. One can reexpress the full 
solution in terms of an initial density matrix p as well: 



p{r,s,t) 



dr'dk 
2nh 



Po(r',S(t)) 



[ ik(r — r ) /•* , , 

x exp I — i— ; - + / di! 

\ n Jo 



iFS(t - Q 
h + T 



G 



's{t-ty 



- 1 



(5.9) 



There are no restrictions on the initial conditions p Q (r,s), but it is convenient in our 
considerations below to use a particular form. If we have an initial Gaussian wavefunction, 
MX) = exp (-X 2 /4a 2 )/(27Rx 2 ) 1/4 , then: 



1 



Po[r,s) 



V2^F 2 

d (s : k) = exp - 



exp | 


' At 2 + s 2 




8(T 2 






2h 2 


8a 2 



(5.10) 
(5.11) 



B. Attractors and Repellors 

It is clear that the time evolution of p(r, s,t) depends on the properties of S(t). Hence 
the flow of the solution p(r, s, t) can be better understood if we examine the stable and 
unstable fixed points of S(t). In order to discuss further the character of this solution it 
is convenient to use a specific form for the correlator G, so for illustrative purposes we use 
G(x) = exp(— x 2 /2). The fixed points are determined from the condition 

which is plotted in Fig. 1 (top). (The analogous result for the case of a periodic correlation 
function G(x) = cos[x] is shown in Fig. 1 (bottom).) This has a maximum value k given 
by 

*° = fe (513) 



The character of the trajectories S(r) determined by solving Eq. ( |5.7| ) depends on whether 
| A; | > ko or \k\ < ko- Since k is not dynamical, the evolution is only along the s direction. 
The flow lines S(r) are shown in Fig. 1 with the arrows for selected values of k. As one can 
see, the part of the curve between — 1 < s/X < 1 is a line of stable fixed points (attractors), 
while for |s/A | > 1, it becomes a line of unstable fixed points (repellors). When \k\ > k the 
r.h.s. of Eq. (|5.7|) maintains a definite sign so S(t) is either a monotonically increasing or 
decreasing function of time, for any given initial condition 5(0) = s. When — k < k < fc one 



14 



can see from Fig. 1 that there are two types of solutions. Furthermore, since ko depends on 
P, as (3 — > 0, k — > and the character of the dynamical evolution depends on temperature. 
In general the trajectory S(t) can be determined through a simple quadratures 



5(r) 



dx 



M 2MX 



x 



n -l 



(5.14) 



C. Momentum Cumulants — Thermalization 

In order to determine the momentum cumulants we must construct the Taylor expansion 
of the function In d(s, k, t) in powers of s at k = 0. One can see from Fig. 1 that along 
the k = line, all trajectories flow to the origin S = in the long time limit. By changing 
the integration variable in Eq. ( |5.5| ) from time to x = s/Xq, using Eq. ( |5.7| ), in the limit 
t — > oo, where all trajectories have the property that S(t) — > we obtain 

1 2FMX 3 r s / x o ( r 2 \ 

In d(s, 0, t) = In do(S(t), 0) + — JZt L dxexp(^) (5.15) 

2MXI r/x dx 




+ h 2 (3 Jo x 

From the power series expansion of the integrands, one can readily read off all the momentum 
cumulants in the limit t — > oo. In this limit the initial condition become irrelevant as 
lado(S(t), 0) — > 0. In particular the first and second cumulants are: 

2MFX 2 F 

«p» = -prin- = ~ (5 - 16) 

((p 2 )) = j = MT (5.17) 

One can see the physical picture emerging here. In the long time limit the quantum 
particle reaches a terminal velocity determined from the first cumulant: 

2FXI F 

(5.18) 



M .iYIi M 7 ' 

with the definition of the friction coefficient 7 identical to that in the Kramers equation, 
Eq. ( |3.2| ). Furthermore, the kinetic energy of the particle equilibrates to the proper thermal 
equilibrium result: 

<£> = M! + &_ = M^k + 1 r 5 191 

2.\/ 2M 2M 2 2' 1 ' ' 

What is more notable, however, is that the momentum distribution has higher than 
second order cumulants, which increase exponentially with the order of the cumulant. In 
the absence of the linear potential [F = 0) only the even order cumulants are nonvanishing, 
and are given by: 
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2d 



nm-ir'^f ^] • (5.20) 
n h (3 \X J 



The presence of F adds only odd cumulants: 



2n-l 



((P 2 - 1 )) = {-l) n -\2n - 3)!!^ ( A j . (5.21) 

All higher than second order cumulants vanish in the strict classical limit h — > 0. These 
cumulants also vanish in the limit Xq — > oo, which shall be interpreted as a weak coupling 
limit to the thermostat (which is the case in statistical physics). When the coupling to the 
"thermostat" is not weak, in the t —>■ oo limit the function d(s, 0, t) (which is the Fourier 
transform of the momentum distribution) is narrower than a Gaussian and which thus leads 
to an equilibrium momentum distribution with longer tails. This is exemplified in Fig. 
2, where we compare the natural logarithm of the momentum characteristic function Rel. 
( |5.15| ) for the case when M = H — (5 — 1, F = with (from narrowest to widest) Xo = 0.1, 
0.5, 1, 2 and X = oo (the Gaussian limit). The presence of a linear potential does not 
modify the absolute value of the characteristic function, only its phase. 

Naively, one would have expected that the coupling to the thermostat is controlled by 
the magnitude of alone. As one can easily convince oneself however, the coupling between 
the two systems is also controlled by the correlation length X . In the limit X — > oo there 
is no energy exchange between the two subsystems, irrespective of the value of the "coupling 
constant" r^. For X = oo the reservoir never responds to the "external agent" and only its 
excitation spectrum acquires GOE fluctuation characteristics if IN-, po and Kq satisfy certain 
well known requirements. 



D. Coordinate Cumulants and Diffusion 



For the coordinate distributions, we compute the Fourier transform d(s = 0, k,t). Since 
s enters as the initial condition, the solutions which characterize the spatial information are 
all trajectories with initial conditions 5(0) = s = and arbitrary k. From Fig. 1, we can see 
that there are three regions to consider: 1) for \k\ < k$ the trajectories have the property that 
S(t — ► oo) approach the attractor exponentially, 2) for \k\ = k , the trajectories approach 
the attractor as an inverse power law and 3) for \k\ > ko, the trajectories diverge linearly 
in time S(t — > oo) — ► — sign(fc) oo. We shall analyse next the behavior of the characteristic 
function d(s = 0, k, t) in these different regimes. 



1. \k\ < kn 



For small k one can linearize Eq. ( |5.7| ) around the origin and solve the simpler equation: 
,S(r)_ k Wr)_k 



dr M IMXl M 



whose solution is 
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S(r) 



7M 



exp(— 7-r) — 



7M' 



In the £ — > 00 limit, by retaining only terms linear in time we obtain 



(5.23) 



iA; 2FX 2 t 
In d{s,k,t)\ s=0 = lnd (S(t),k) + — 



1 / A-N z AV^X 2 t 



(5.24) 



For an initial gaussian wavepacket, \nd (S(t), k) = — S 2 (t) / &<r 2 — k 2 a 2 / 2fi 2 . In the long time 
limit, S(t) approaches the fixed point Si(k) of (5.12), where Si(k) = —kM/'y + o(k 3 ). From 
this we can determine the first two spatial cumulants: 



((r 2 )) 



2FX 2 



a 2 + 



ti 2 



4M 2 a 2 -f 2 (3 2 Tth 



AX 2 9 
f t = r 2 + 2Dt. 



(5.25) 
(5.26) 



The constant term in ((r 2 )) has two contributions. The first term is due to the initial width 
of the gaussian, while the second emerges in the long time limit, and is absent at £ = 0. 
Physically we see a consistent picture of the quantum dynamics. The particle position grows 
linearly, with the velocity given by the terminal velocity obtained from ((p)). Furthermore, 
the average position displays diffusion consistent with Brownian motion, which can be used 
to determine the diffusion constant D. This is the same expression for D as obtained from 
the fluctuation-dissipation theorem in the classical limit (Kramers equation) in Section III. 

As with the momentum distribution, the coordinate distribution in not Gaussian, and 
has longer tails. (In the Brownian motion limit, these tails vanish; see below.) An analytical 
explicit construction of the entire spatial distribution and its cumulant expansion is not 
quite trivial, as the characteristic function d(k, 0, £) has singularities and different asymptotic 
time behavior depending on the value of k. In particular, for small values of k the position 
of the repellor 82(h) 



sign(A;)Xoy — 2 In \k\ is not analytical around k = 0. This is an 
indication that either various moments of the spatial distribution do not exist (perhaps they 
are divergent, as in the case of the Cauchy distribution) or they increase with time at a rate 
much faster than linear, similar to the case (3 = discussed in Ref. ||. In other words, the 
function In d(s, k, t)\ s= o has some singularities in the A;-plane, which have to be dealt with 
more carefully. 

More generally, for \k\ < k and large times Eq. (|5.7|) can be written approximately as 



dS(r) 
dr 



2MXl 



exp 



si(kf 
2X 2 



si(kf 
X 2 



(S(r)- Sl (k)), 



(5.27) 



where S\(k) is the A;-dependent position of the attractor and the boundary condition is 
in this case S(t — ► 00) — > s\(k). In computing the leading term in t one can use this 
approximate trajectory. In this way one arrives at 



lnrf(s, k, t)\ 



s=0 



In d (s 1 (k), k)\ s=Q + 



iF Sl (k) T l 



exp 



' 2X 2 



- 1 



(5.28 



17 



For small values of k the function d(0, k, t) is narrower than a Gaussian, as one can establish 
easily by comparing the Taylor series in k of this expression with Eq. ( 5.24Q . Since the 
Fourier transform of this function is nothing else but the spatial distribution, we can thus 
conclude that at large distances the spatial distribution has longer tails than a Brownian 
particle. 



2. 



kn 



For this critical value of k, the trajectory with initial condition s = 0, k = ±/c , will 
approach S = =FXo in the limit r — ► oo, as seen in Fig. 1. It is sufficient to look at the 
case k = ko- To examine the behavior of S(t) in the neighborhood of the fixed point, it is 
convenient to take: 



S(t) = -X (1-s(t)). 



+ 7(1 - e) exp 



The dynamics is then given by 

cfe(r) k 
dr ~ ~MX 

The solution in terms of S then has the power-law behavior: 

k r 



MX 



S(t) 



-Xn 



k r + MX 

Substituting this into the solution for d(s = 0, k, t), we obtain: 



lnd(0,fc ,t) = lnd (S{t),k ) + 



k 



MX 7 



1 h 



FX 



+ 



FMXn 

h I - 

h'y fik 



In 



kn 

1 + — —t 

MXn 



+ o(k, 



(5.29) 



(5.30) 



(5.31) 



(5.32) 



3. \k\ > k 

In the other regime, when \k\ 3> k and t — ► 00, j3 no longer plays a role. Strictly 
speaking, this analysis is only valid in the (3 = limit. The approximate equation for the 
trajectory is 



dS(r) 
dr 



k 

M' 



(5.33) 



and thus the characteristic function acquires approximately the form 



lnd(s, k, t)\ 



s=0 



\nd (—*)- 



iFkt 2 T l 
2hM + T 



2nMX t 



2\k\ 



-t 



(5.34) 
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Notice, that if for \k\ < ko the characteristic function is narrower than a Gaussian, for 
| A; | > k one has just the opposite, the characteristic function tends to a time-dependent 
constant (modulo a nontrivial phase however). The above relation shows that the short 
distance behavior of the propagator has a quite unusual behavior. In order to construct the 
propagator one should use the initial density distribution po(r, s) = S(r) and thus do(s, k) = 
1). The above behaviour for large k implies that even at finite times an exponentially small 
part of the initial spatial distribution is left at the origin, namely exp(— FH/h)S(r). 



4- General Behavior for all k 

An alternate manner to construct the coordinate distribution d(0, k, t) is by solving 
numerically Eq. (|5.7| ) for S(t) and substituting the trajectory directly into Rel. (|5.5| ) for 
d(0, k, t). In Fig. 3 (top) we compare the numerical solution to Eqs. (|5.5| ) and ( |5.7|) 
(solid) with the Gaussian limit obtained by keeping only the second cumulant (dashes). The 
solutions are obtained for a time t = 100, and Xo = 1. The results are all similar. However, if 
one goes to shorter correlation lengths, the importance of higher order cumulants is striking. 
In the bottom figure we show the same for a much shorter correlation length, Xo = 0.1, also 
at t — 100. All the remaining parameters are kept at unity. As one can see, the exact result 
drops off abruptly at k — k . The reason is that for \k\ > k , the function S(t) escapes 
towards — 00 linearly in time while the solution for \k\ < ko slowly converges to the stable 
fixed points. If one factors out the overall exp(— TH/Ti) from \d(s, k,t)\, the reminder grows 
exponentially in time for \k\ < ko and tends to a time independent function for \k\ > ko and 
as a result the discontinuity at \k\ — ko becomes thus more pronounced. 



E. Brownian Motion Limit 



Even though we have shown that the system equilibrates to the correct thermal limit, 
the time evolution towards this equlibrium state is rather complex. For k = all trajectories 
have the same asymptotic behaviour, S(t — ► 00) — > 0, irrespective of the initial conditions. 
It is not difficult however to see that if a trajectory starts far away from the origin, it will 
take an exponentially long time to reach the neighborhood of the origin, as the r.h.s. of Eq. 
(d77) is exponentially small for \s\ ^> X . 

One can see from the expressions for the cumulants that if the limits (see also Ref. ||) 

0, (5.35) 
(5.36) 



h 

Xo 
Xo 



are taken, with the friction coefficient 7 remaining finite, one obtains the case of pure classical 
Brownian motion. All but the first two cumulants for coordinate and momenta vanish, and 
one is left with a Gaussian process. These limits can be achieved also by keeping H finite 
and thus obtaining the case of a quantum Brownian particle. 
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VI. QUADRATIC POTENTIAL 



A classical particle in a harmonic oscillator potential, treated with the Fokker-Planck 
or Langevin equations, will thermalize, with the equilibration given by the virial theorem. 
When the particle is quantum, the fluctuations are chaotic and moreover the coupling be- 
tween the two subsystems is finite, we observe significant departures from the this idealized 
situation. 



A. Exact Solution 



For a particle in a harmonic oscillator potential, we start with 

H Q (X) = 



Ti 2 o2 Mlu 2 X 2 



(6.1) 



2M A 2 

The equation for d(s, k, t) can be solved through quadratures, in the same manner as for 
the linear potential. We note that the solution we obtain is not limited to the simple case 
under consideration. Analytic solutions are also possible if we want to include a linear 
time independent potential term, a general quadratic time dependent potential, and/or a 
multi-dimensional treatment. 

We shall look for a solution using the same representation of the density matrix intro- 
duced in the previous Section in Rel. ( [2.17 ). The equation for density matrix is in this 
case 



inot + —o r o s , , 

M 2X M 



Xq 



p(r,s,t) 



(6.2) 



, ;r l 



G 



s 

X~ 



1 



p(r,s,t) 



and the corresponding equation for the function d(s, t, k) is: 



d t + 



M 2MX \X 



d s - Muo 2 sd k \ d(s,k,t) 



h 



G 



s 

X~ 



d(s,k,t). (6.3) 



The solution of this equation can be obtained again using the method of characteristics [15 



In this case we will have a two-parameter solution for d(s, k, t) which will depend on the 
functions S(t) and K(t) which satisfy the auxiliary equations: 

(3T l h 



dS{r) 

dr 
dK(r) 

dr 



K{t) 
M 



Mlv 2 S{t) 



2MXn 



-G' 



S(r) 



with the initial conditions 



S(t = 0) = s, and K(r = 0) = k. 



(6.4) 
(6.5) 

(6.6) 



We do not write the full solution here since it is more convenient to solve the equations in 
action-angle variables. Before we do so, let us examine the fixed point structure which will 
emerge. 
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B. Attractors and Repellors 



The flow in the "phase space" (S, K) is not Hamiltonian in character, which is not 
surprising. Some typical trajectories are shown in Fig. 4 for G(x) = exp(—x 2 /2). In this 
case there is no qualitative difference between the dynamics of this correlator or a periodic 
function such as G(x) = cos(a;). In the top of Fig. 4 we have lo — 1, 7 = 0.1 while in the 
bottom the friction is stronger, with lo — 1, 7 = 0.5. The general pattern of the trajectories 
in the (S, K)-plane seems to be quite simple, in the limit r — > 00, irrespective of the initial 
conditions, all trajectories spiral counterclockwise around the origin (unless the motion is 
overdamped). The origin is thus a stable focus. There are no trajectories going away to 
infinity in any direction in the "phase space" (S,K). The plane (S,K) is separated into 
four regions by the two curves 



S = and 



K 
M 



2MX 



a 



s_ 



0. 



(6.7) 



Each of these lines correspond to K(t) = and S(t) = respectively. Near the focus it is 
simpler to solve the linearized equations of motion 



dS(r) 
dr 

dK(r) 
dr 



K(t) 



T 



M 

= Muj 2 S(t), 



IMXl ' 



which can be solved analytically. If the condition 

(3T l h 7 



4MX 2 

is fulfilled, then the required trajectories are 



(6.8) 
(6.9) 

(6.10) 



S(t) 
K(t) 



2k + M'js 

s cos cur — sin cur exp 

2MW 

-fk + 2Mu 2 s . 

k cos lot H — sin lot 

2lo 



7T\ 

' 2 ) 



7T 

exp l - — 



(6.11) 
(6.12) 



where 



LO = 



( pnn 



4MX 2 / 



= I w 



T' 



(6.13) 



The case 00 < 7/2 is formally similar to the case lo = iQ to be discussed in the next section. 
Strictly speaking, in the first order in j3, in which we have derived our formulas so far, one 
has lo ~ lo and we shall use mostly this approximation from here on. 



21 



C. Solution in Action Angle Coordinates 



The flows shown in Fig. 4 suggest that action-angle coordinates might be better suited 
to the dynamical analysis of this problem. The analytical construction of the higher order 
cumulants is somewhat cumbersome and we shall limit our analysis to some general features. 
The action-angle variables are in this case 



2Mu 



arctan ■ 



Mus 2 



k 



Mus 



(6.14) 
(6.15) 



Here / is related to the energy. We can now rewrite the evolution equation Eq. (|6.3|) in 
these variables 



d t - ud^ + 7 exp 



I cos 2 



sin cos 



ri 
h 



exp 



+ 21 cos 2 
/cos 2 



MuXl 

The equations for the trajectories acquire the following form 




(6.16) 



d$(r) 
dr 

dl(r) 
dr 



7 

u — — sin 2<3>(t) exp 



-7l(r)[l + cos2$(r)] exp 



2l^ 1 + COS2 ^ 
2MluX 2 1 



:[1 + cos2$(r)] 



(6.17) 
(6.18) 



with initial conditions 



1(0) 
$(0) 



Mus 2 



2Mu 



arctan ■ 



Mus 



(6.19) 
(6.20) 



Irrespective of the initial conditions the action I(r) is always a monotonically decreasing 
function, vanishing in the long time limit. If the motion is not overdamped (i.e. u > 7/2) the 
phase $(t) is a monotonically increasing function of time. In the weak friction limit, when 
7 < w, one can replace the equations for the trajectories with the following approximate 
equations 



d$(r) 
dr 

tfl(r) 
dr 



u, 



-7I(t) exp 



Kr) 



2MuX 2 



Zn 



T 



2MUXZ 



Xi 



2MuX 2 0/ 



(6.21) 
(6.22) 



obtained after averaging the initial equations over the fast motion (i.e. over one period 
2-k/u). Here Tq{x) and 1\{x) are the modified Bessel functions of first kind. The r.h.s. of 
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Eq. ( |6.22j ) behaves linearly in I for small values of the action and as 1/ y/l for large values. 
This behaviour is consistent with the fact that the friction term is effective only near the 
origin in coordinate, i.e. for s ~ O(X ). 

After manipulations similar to those used in the previous section, and averaging over the 
fast motion and by changing the integration variable from time to action (using Eq. ( |6.22| )), 
one obtains that the asymptotic expression for the characteristic function is 



\nd(I,(f),t) 



2MXl 
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2MuXl 
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2MuX^ 



(6.23) 



2MwX 2 



where / = 1(0), which was defined above as the initial value of the action variable. One 
can construct in a straightforward manner all the spatial and momentum cumulants, by 
reverting to the initial space-momentum variables s and k. Note that d(s, k, t) = d(I, <f), t) 
is time independent in this limit, as expected. One can extract easily the behaviour for large 
and small /. 



M 2 



\nd(I,(f),t) 



2ixMX 2 _ 



3/3h 2 



3/2 



I > MuXl 



(6.24) 



By taking the Fourier transform of the above expression for d(s,k,t) for either k = or 
s = one can determine either the momentum or the spatial equlibrium distribution of a 
harmonic oscillator coupled with a "reservoir". For the oscillator, the action variable I is, 
up to a trivial factor, simply the total energy E. Hence one can also extract the energy 
distribution from the above expressions. 

D. Eigenvalues of the Time Evolution Operator 



It is instructive to construct the eigenvalues and the eigenvectors of the time evolution 
operator. If we rewrite the time evolution equation for the density matrix in the form, 



ihd t p(r, s, t) = Op(r, s, t) 
we can consider the eigenvalue problem associated with O: 

Op(r, s) = Ap(r, s). 



(6.25) 



(6.26) 



The equilibrium solution will correspond to A = 0, and in general the spectrum should be 
complex. While an analytical solution does not seem possible for a general G(x), if we take 
G(x) = 1 — x 2 /2, one can readily solve the problem. 
For the equilibrium state (A = 0) one obtains 



Po(r,s) = exp 



(3Mu 2 r 



2^2 



M 



(6.27) 
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We will discuss the physical properties of this solution in the next section. This problem 
now becomes identical to the eigensolutions of the Fokker-Planck equation for the oscillator 



TM. Following that analysis, the eigenvalue spectrum is given in terms of two integers, ni 



and n 2 . The basic roots are 



\ ± = -^±nuJ, uJ=]Juj 2 --^ (6.28) 

and the entire spectrum can be written as: 

V,n 2 = — — (ni + n 2 ) + ^(rii-n 2 ), (^ = 0,1,...) (6.29) 

The time evolution of the density matrix will equilibrate to the A = eigenvector, as all 
other components will decay in time as exp[— (n\ + n2)'~ft/2h]. The temperature dependence 
appears indirectly through 7, which vanishes in the f3 = limit. These patterns are shown in 
Fig. 5 for (3 > 0, where it is clear that all eigenvalues with the exception of the equilibrium 
state (A = 0) lead to decay. For (3 = 0, the spectrum collapses to the real axis, and all points 
A = nuJ are infinitely degenerate. 



E. Recovery of Equilibrium Thermodynamics 



For the harmonic oscillator, we would like to see whether or not the random matrix bath 
can act as an ideal heat bath. The quantum equilibrium density matrix we would expect is: 



Peq(X, Y) = ^ exp 



n=0 



-P[n+~)?uj 



<P n {X)<p n {Y) 



where n are the oscillator wavefunctions given by: 



<Pn{X) 



Of 



ir2 2n (n\) 2 



1/4 



a 2 X 2 * 



H n (aX) exp 



a 



Muj 

IT' 



and H n is the Hermite polynomial. Defining z = exp[—/3hu], we can write: 



Pe q (X,Y) 




exp 



a 2 (X 2 + Y 2 



az 



tt(1 -z 2 ) 



exp 



a 



2(1 - z 2 ) 



00 ( z\ n 1 

E ff) -H n (aX)H n (aY) 

J n=0 VZ/ U - 



(X 2 + Y 2 )(1 + z 2 ) -AXYz 



(6.30) 

(6.31) 

(6.32) 
(6.33) 



In the last line we have used the generating function for Hermite polynomials. If we take 
the leading order in {3 and transform from X, Y to r, s, we find 

> a 2 (3Tiuo 



Pe q (r,s) 



a: 



2n(3h 



exp 



Q- 



2f3htu' 



(6.34) 



which is precisely the A = eigenvector found in the previous section. Hence, to leading 
order in (3 (and % in a sense, as (3 enters as (3%), and in the Brownian Motion limit, we 
recover the equilibrium density matrix. In general the equilibrium density matrix will have 
a non-Gaussian character. 
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F. Cumulant Expansion 



To get the cumulants, one can use the procedure described previously (see also Ref. |§) 
and derive that in the t —>■ oo limit 

«p 2 » _ ((M,V)) _ 1 _ T 

2M ~ 2 _ 2(3 ~ 2 ' [b - 6b) 

as one might have expected. The higher order cumulants can be obtained by setting k = 
or s = respectively in the Taylor expansion in I of the Rel. ( 16.231) . We shall not try 
here to derive explicit expressions for higher order cumulants, as they apparently seem to be 



quite complicated and hardly very revealing. As Rel. ( 3.24 ) clearly shows, the equilibrium 
distribution is much narrower than a Gaussian one, which is already indicative of the presence 
of nonnegligible higher then second order cumulants. In Fig. 6 we compare lnd(I, <f),t —* oo) 



as a function of \fl obtained from Eq. ( |S.23[ ) for X = 0.5, l,2,oo (narrowest to widest). 



The widest X = oo is the equilibrium distribution exp(— [3H ), (the first term in the low 
/ expansion, see Rel. ( |6.24 )). Notice that the equilibrium distribution is independent of 



<p as well as of the actual values of the friction coefficient 7 and of the coupling constant 
T^. Except for a "trivial" overall factor, the shape of this function is controlled by a single 
parameter, the "characteristic action" 2MujXq, which depends on the "roughness" of the 
coupling to the "reservoir" . 



VII. INVERTED PARABOLIC POTENTIAL 

The tunneling through a potential having the shape of an inverted parabola can be 
studied by performing the formal replacement u = iflo in the corresponding equations of 
the previous section. 



A. Fixed Points 

The equations for the trajectories (S(t), K(t)) become now 
dS(r) 



dr 
dK(r) 



K(r) (3T±h Ql 




(7.1) 



dr 

The two curves we have introduced in the previous section 



M 2MX 

MQ 2 S(t). (7.2) 



S = Q an d £-J^L^(4r)=0, (7.3) 



M 2MX \X 

play a similar role in this well. However, the general flow pattern in the "phase 

space" (S, K) has an entirely different aspect, as shown in Fig. 7. The origin is now an 
unstable focus point. Near the origin we can again use the linear approximation G'{x) ~ — x 
and thus consider the much simpler equations 



25 



dS(r) _K{t) _(3T^hS(r) 

dr ~ M 
dK(r) 



2MXl 



dr 

The imaginary eigenfrequencies are 



-MQ 2 S(t). 



(7.4) 
(7.5) 




(7.6) 



with Q + > and f2_ < 0. One can easily recognize that fl + is the renormalized attempt 
frequency in the classical Kramers result for the transmission through a barrier 0]. 

As we have noted in the previous section the case of an inverted parabolic potential 
bears some formal similarities with the case of over damped motion (u < 7/2), with only 
one difference, that both imaginary eigenfrequencies are negative. We shall not discussed 
this case in detail in this paper. 

For arbitrary initial conditions S(0) = s, K(0) = k in the neighborhood of the origin 
almost all trajectories are diverging exponentially away as exp(fi + r), except when s and k 
are exactly along the line 



1 - — ) MQs. 
20/ 



For sufficiently large s the stable manifold is defined by the line 

k = M£l s. 



(7.7) 



(7.8) 



For small friction coefficient (7 <C Qq) these two lines are almost identical. For large values 
of the friction coefficient however (7 3> f2o) the stable manifold near the origin almost 
coincides with the Os-axis. 

The solution to the above linearized equations of motion are 



S(r) 



1 - 



+ 



K{t) 



+ 



=: 1 + 



1 + 



1 - 



1\ 


- 


k 


2QJ 




2MVL 


1_\ 


+ 


k 


2QJ 


2MVL 






Mn 2 s 


2QJ 


- 


2VL 


7 > 


+ 


Mn 2 s 


2QJ 


2VL 



exp 



exp 



exp 



exp 



n 



(7.9) 



n 



+ l lT 



(7.10) 



With increasing time and for initial conditions such that IS'(O) = s\ ^> X , the term re- 
sponsible for friction in Eqs. (|7.2j ) becomes negligible and the trajectories are governed by 
exponentials exp(±f2 T). 

In the case of very small friction (7 Q ) one can obtain an approximate expression 
for the trajectory S(t). If one is interested in the character of the spatial distribution and 
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the way this is modified in the presence of dissipation, one needs trajectories with initial 
conditions 5(0) = 0, K(0) = k only, namely 



S(r) 
K(t) 



k sinh UnT 



Mtt 
k cosh Qqt 



cxp 



o 2 



k^{r) sinh Q r 



exp 



[ T dr' 
Jo 



7 (r' 



where the "instantaneous" friction coefficient is defined as follows 

k 2 sinh 2 Q r\ 



7(t) =7 exp 



2M 2 fi 2 1 X 2 



(7.11) 
(7.12) 

(7.13) 



Using these expressions one obtains for the characteristic function for the spatial distribution 
the following expression 



d(0, k, t) = d (S(t), K{t)) exp | jf dr 



cxp 




- 1 



(7.14) 



In spite of the fact that it is hopeless to evaluate analytically this expression as a function 
of k and time t (there are at least five exponentials nested into one another) , even a cursory 
analysis of the /c-dependence originating from dissipation alone shows that this function is 
significantly narrower than a Gaussian. This at once is indicative of the presence of higher 
than second order cumulants in the spatial distribution, but most of all of the fact that the 
spatial distribution has a wider spread than the corresponding one obtained in traditional 
approaches. If one were to perform a Taylor expansion of the integrand in powers of k to 
determine the cumulants for example, one observes that each power of k is always multiplied 
by sinhf2 T. Upon performing the time integral one thus obtains that all cumulants increase 
exponentially with time and that in particular all of them are positive as well. 



B. Eigenvalues 



The eigenvalues of the evolution equation are obtained from those of the harmonic oscilla- 
tor, by making the transformation uj = ifl . Again, specializing to the case G(x) = 1 — x 2 /2, 
we obtain the spectrum: 



A+ = — %h 



± Jng + 



r 



(7.15) 



A 



i'yh 



9 v m + n 2 ) +ihJSl% + ^-(rai -n 2 ), 



ni,n 2 



n. ; 



1,2,...) 



In this case the eigenvalue A = 0, which corresponds to an equilibrium solution, is no longer 
present. The eigenvalues are all purely imaginary, and have values above and below the real 
axis. The occurrence of eigenvalues with positive imaginary parts is due to the unphysical 
asymptotic behaviour of an inverted parabola potential, and thus to the absence of an 
equilibrium distribution. 
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C. Remarks on Tunneling 



Kramers showed that dissipation leads to a reduction of the flux through an inverted 
parabolic barrier, as the unperturbed attempt frequency in the transition state theory should 
be replaced with the renormalized one Qq —>■ Q + < Qq One basic assumption in Kramers 
approach was the fact that the friction coefficient is momentum independent. The present 
approach can be interpreted as a theory with a momentum dependent friction coefficient, 
which in the zero velocity limit reduces to the classical value. For finite velocities however, 
the effective friction coefficient in our approach is smaller than the one for zero velocity (see 
the above approximate solution for the trajectory for the case of small friction). One can 
thus expect two effects: i) the effective or average attempt frequency in our approach is in 
between the Kramers value and the value corresponding to no friction, i.e. fl + < Q e ff < ^o! 
ii) the spatial density distribution is also modified. Unfortunately the rest of our argument 
concerning the behaviour of the spatial distribution is pure qualitative at the present time. 
In the classical transition theory the flux through the barrier is controlled by the probability 
population of the states at the top of the barrier. With decreasing friction coefficient, the 
thermalization of the spatial distribution in the classically allowed region is less effective and 
as a result the overall flux through the barrier decreases as well. In the previous section we 
have established that the energy distribution has generally wider tails than a pure Boltzmann 
distribution and thus states with higher energies have a higher relative probability to be 
populated when the coupling strength between the "reservoir" and the subsystem is finite. 
It is thus reasonable to expect that tunneling states would have a higher relative occupation 
probability in our approach. Overall the effect of an effective momentum dependent friction 
coefficient is likely to lead to an enhancement of the tunneling probability when compared 
with the classical Kramers result. Since the role of the pure quantum effects and the role 
of the finite coupling strength to the "reservoir" are intimately related, we cannot at the 
present time disentangle the specific role played by each one of them separately. 



VIII. TUNNELING IN A DOUBLE WELL POTENTIAL 

For the problem of tunneling, we consider also the dynamical evolution of a particle in 
the double well potential given by: 

H (X) = -— dl + a(x 2 -— ) . (8.1) 
ov ; 2M x \ 2a J K ' 

For this or more complex potentials, analytic solutions are no longer possible. The form of 
the potential is shown in Fig. 8, together with the first few eigenvalues corresponding to the 
parameters a = 5, b = M = 1/2. 



A. Tunneling Rates 

The numerical solution of the evolution Eq. ( |2.12p for the density matrix is done as 
explained in Sect. IV.B. Let us examine first the effect of dissipation on tunneling. For this 
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purpose, we first put (3 — and solve the equation for some representative values of the 
parameter 1^. Let (f>i(X) and 4>2{X) be the first two eigenstates of the double well. We take 
as initial state the combination 



(8.2) 



which represents a wave packet that is mostly localized in the left well of Fig. 8. In view 
of the fact that splitting between the first two states is very small, AE ~ 3 x 10 -6 , the 
tunneling of the wave packet to the other well would take a very long time should there be 
no dissipation. Fig. 9 shows the strength that tunnels through the barrier as function of 
the time for some values of TK At the beginning of the tunneling process, the rate increases 
almost linearly with this parameter. For very large time, all the tunneled strengths approach 
the limiting value 1/2 corresponding to the equal distribution between the two wells. 



Next, we study the dependence of the tunneling on the temperature. For this, we fix a 
value of r-L and solve the evolution equation for different values of (3. One would expect that 
the rate would increase with the temperature, i.e. when f3 decreases. Fig. 10 shows such a 
behaviour for five values of f3 — 0., 0.05, 0.10, 0.15, 0.20. For larger values of (3, the strength 
on the right hand side may come out negative for reasons we discussed in Sect. II. 



At this point, we are also able to compare our results with that obtained from the 
Caldeira-Leggett model. To reproduce the latter, it is sufficient to replace the Gaussian 
form of the correlation function that we have used up to now by its quadratic approximation 
G(x) = 1 —x 2 /2 .The effect may be seen better by choosing a smaller value for the correlation 
length X . Fig. 11 shows the density distribution after the time t — 4 for the parameters 
X = 0.2, (3 — 0, — 1. The difference is quite appreciable (notice the ordinate logscale), 
not only in the tunneled strengths but also in their shapes. 



We have numerically computed the eigenvalue spectrum by expanding the density matrix 
in the eigenvectors of this quartic potential. Taking N max = 31, we have diagonalized the 
evolution operator O. The complex eigenvalues with the smallest imaginary parts are shown 
in Fig. 12 (top) for (3 = and 0.2. If we denote the eigenvectors as pt(X, Y), then the time 
evolution, which includes the tunneling, has the form 



B. Temperature Dependence 



C. Comparison to the Caldeira Leggett Approach 



D. Eigenvalues of the Evolution Operator 




(8.3) 
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One can see in Fig. 12 that the tunneling can be influenced by several nearby eigenvalues, 
and that as (3 increases, the patterns of important eigenvalues changes considerably. In 
the bottom of Fig. 12, we plot the imaginary parts of the eigenvalues as a function of 
temperature for the ones with the smallest imaginary component. One observes that the 
imaginary part of the eigenvalue with the smallest nonvanishing absolute value apparently 
vanishes in the limit T — > 0. This is consistent with Caldeira and Leggett's conclusion 
that at T = dissipation decreases the tunnelling rate. This is not true however at high 
temperatures, as one can see. 



IX. CONCLUDING REMARKS 

We have developed a dynamical theory of simple quantum systems coupled to complex 
quantum environments, where the environment is a general "chaotic" bath of intrinsic exci- 
tations. The model Hamiltonian we introduce for the intrinsic subsystem incorporates the 
generic properties of finite many-body systems. This includes an average level density of 
states sharply increasing with energy, the presence of universal spectral (or random matrix) 
fluctuations for the intrinsic system and the variation of these properties while changing the 
"shape" of the intrinsic system modeled with parametric banded random matrices. In this 
way, the intrinsic system is capable to easily absorb energy due to its large heat capacity. 
We did not yet allow the intrinsic system to perform mechanical work, but this feature can 
be readily implemented. 

We are able to solve the dynamical evolution equations for these systems, without mak- 
ing any uncontrollable approximations or assumptions. Using the Feynman- Vernon path 
integral approach, we derive the random matrix influence functional, and an evolution equa- 
tion for the density matrix of the "slow" subsystem at finite, but high temperatures. This 
evolution equation is surprisingly easy to manipulate and in many instances one can con- 
struct full solutions by quadratures, for cases when the corresponding path integral can be 
computed only by brute force. The entire treatment is performed at the quantum level. We 
have shown on the other hand that in the classical limit, the evolution equation reduces to 
the Kramers form. 

Our analysis was limited to the case when the motion of the simple system can be treated 
in the adiabatic approximation. Thus, we have not taken advantage of another parameter 
in our description of the reservoir, the bandwidth of the random matrix k . It is well 
known that banded random matrices lead to localization |18||. It would be very interesting 



to explore whether localization in the "reservoir" states induces localization in the simple 
quantum system coupled to such a reservoir. (One can show that for finite bandwidths k 
the influence functional in our model can be brought to a form similar to the one suggested 



by Cohen |TT| to describe localization, although we did not analyze this limit here.) 

The quantum evolution equation we have derived for the density matrix is not more 
complicated than a classical Fokker-Planck equation. However, this evolution equation 
describes processes which first of all are quantum in nature and moreover have manifestly 
non-Gaussian features. So although we started with a Gaussian process characterized by 
our parametric random matrices, we ended up with a non-Gaussian dynamical process. 
Unlike traditional approaches, we have not assumed that the coupling to the "reservoir" is 
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vanishingly small. In particular, this is one of the reasons why the diffusion process acquires 
a non-Gaussian character. The spatial and momentum distributions are characterized by 
high order cumulants, which increase exponentially with the order. As we have discussed 
earlier ||, this is indicative of the fact that the corresponding distributions have longer 
tails than previously expected. In order to put in evidence such features, one ought to be 
careful about the order in which the various limits are taken. If one takes the semiclassical 
limit first in the evolution equation for the density matrix, such features are lost. If in the 
actual quantum solution one takes either the semiclassical limit (h — ► 0) or the limit of weak 
coupling to the reservoir (either IN- — > or/and X — ► oo) then one recovers the standard 
results as well. We have thus shown that, at least in the framework of the present model 
(which we believe to be sufficiently generic), both semiclassical and weak coupling limits 
have a rather singular character, which is seen in the temporal domain or/and the phase 
space distributions. 

The particular form of the quantum evolution equation we have discussed in this paper 
has one rather unneccessary restriction, based on the "translational invariant" form of the 
influence functional we have used, see Eq. ( [2. Ill ), which lead to Eq. ( |2.12| ). One can easily 



convince oneself that this form we have assumed can be easily amended and, to within some 
approximation, a more general quantum evolution equation for the density matrix of the 
"slow" subsystem can be derived, namely: 

ihd tP (X,Y,t) = {^~^ + U ( X ) - U ( Y ) (9- 1 ) 



7XjdG(X,Y) 
2 ds 



[P x - Py] + iT* [G{X, Y) - 1] j p(X, Y, 



t) 



where s = X — Y and where the correlator G(X, Y) satisfies the only symmetry require- 
ments G(X,Y) = G{Y,X) = G*(X,Y). One has to interpret 7, /3, T l and G(X,Y) as 
phenomenological quantities, characterizing the "reservoir" and its coupling to the "slow" 
subsystem. We thus have now at our disposal a quantum evolution equation, with an effec- 
tive velocity dependent friction coefficient, a finite coupling strength to the "reservoir" and 
which in the classical and weak coupling limits reduces to the Kramers equation. In the 
weak coupling limit only this evolution equation is equivalent to the high temperature limit 
of the Caldeira-Leggett model. 
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APPENDIX A: EVOLUTION EQUATIONS FOR THE AVERAGE 
PROPAGATOR AND THE INFLUENCE FUNCTIONAL 



The time evolution of the fast subsystem is found by solving the time-dependent 
Schrodinger equation in the form USHlQfl: 



<f>(t) = Texp 



<Ko)=w(x(t))0(o). 



(Al) 



where T is the time ordering operator, and U(X(t)) the propagator. (We assume that the 
initial state 0(0) is uncorrelated with the Hamiltonian Hi(X(t)) at later times; correlated 
initial conditions have been discussed elsewhere 0.) One can show that in the leading 
order in an expansion in 1/Nq the average propagator U(X(t)) = U(X(t)) is diagonal in the 
representation we have chosen. Its diagonal matrix elements have the following form 



U k (X(t)) = (k\Tex V 



~ I dsHi(X(s)) 
n jo 



|k) =^(--f-)cr(X(i)) 



(A2) 



(note that a(X(t)) is state independent) and a(X{t)) satisfies the following integral equation 



r\ r s \ , , , fx(si) -x(s 2 ) 

a(t, t ) = 1 - — / dsi / ds 2 a{s 1 , s 2 )a{s 2 , t ) "/H^i - s 2 )G 

tl Jto Jig \ v\q 



(A3) 



In the following we shall consider that to = 0. Here Pp(s) is the fourier transform of the 
matrix band form-factor in the correlator [/i 1 (X)]j :; [/i 1 (y)]^, defined in Eq. 



Pais) = PS(-s) 



2nh 



exp 



K 



2h 



°- 2 [s + M 
2 \ 2 



(A4) 



The influence functional can be determined by solving the following evolution equation: 

C(X(h), Y(t 2 )) = a(h, 0)a*(t 2 , 0) + L ds x t ds 2 C(X( Sl ), X(s 2 )) (A5) 

n Jo Jo 



xP*J Sl -s 2 )G 



, X(s 1 )-Y( 82 ) 



X n 



(T(t 1 ,si)a*(t 2 ,s 2 ). 



APPENDIX B: AVERAGE PROPAGATOR AT FINITE TEMPERATURES 



For a given path X(t) the averaged propagator a satisfies the following equation 



<r(Mo) = 1 



dsi / ds 2 a{s 1 ,s 2 )a(s 2 ,t )Pp(s 1 - s 2 )G | - 



h 



to J to 



x n 



(Bl) 



where 
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Pft(s) 



K 



2-kTi 



exp 



2K 1 



s + % 



hp" 



2 as 



(B2) 



In the above equations o"(s!,s 2 ) represents the averaged propagator from time s 2 to time 
Si, as for an arbitrary path X(t) it is not obvious that this propagator depends only on the 
time difference si — s 2 , as was the case in the adiabatic limit and for (5 — 0. We shall assume 
that the following replacement 



Po(s) 



K 



2txU 



exp 



2h 2 



6(8) 



(B3) 



is legitimate and compute the first correction in (3 to the propagator. 

In the equations for the propagator and the influence functional Pp{s) enters under 
integrals with some arbitrary functions as follows 



/ dsP p {s)F{s) » / dsP (s)F(s) + ^ / 



iph [0 dP (s) 



+ — I ds """ K " J F(s) 
2 J-oo ds K ' 



1f ( o) + ^Lf(o) ^ d JM 

2 v 1 2V2^ A ds 



(B4) 



s=0 



where we have used the following relations 

f dsP (s)F(s) = ±F(0), 



J —<x ds 



~~ m ~ 2 — 



Po(0) = 



2nh 

K 



s=0 



2nh 



After a few simple manipulations and after taking into account that 



a(s,s) 



da(s 1 ,s 2 ) 



dsi 

d_ Q f x( Sl )-x( S2 y 

ds 2 V 



da(s!, s 2 ) 



Sl=S'2 



dSn 



0, 



Sl=S 2 



= _± G (X(si)-X(s 2 y 



dsi 



= o, 



(B5) 
(B6) 
(B7) 

(B8) 
(B9) 

(BIO) 



one can easily show that up to terms of order 0(f3 2 ) the averaged propagator is given by 
the following expression 



a(t, 0) = a(t) = exp 



~2h [ + ^~^~ 



(Bll) 



In the equation for the influence functional we shall need the quantity 
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a(t)a*(t') = exp 



exp 



T^(t + f) f0Ko i/3T l \T l (t-t') 



2h 
r^t + t') 



+ 



^7T 4 



2ft 



2ft 



(B12) 



We are allowed to neglect the imaginary contribution in the exponent as in all relevant 
integrals t -t' = 0(J3). 



APPENDIX C: INFLUENCE FUNCTIONAL AT FINITE TEMPERATURES 

In order to simplify somewhat the derivation we shall introduce the "reduced" influence 
functional at temperature (3 



A^(X(ti),y(* 2 )) = £ /3 (X(t 1 ),y(t 2 ))exp 



r+(*l + t 2 ) 

2ft 



(CI) 



which, up to correction terms of order 0{(3 2 ) in the averaged propagators, satisfies the 
equation 



A (X{t 1 ),Y(t 2 )) = l 



(C2) 



^ jy Sl jys 2 Ap(x(s 1 ),Y( S2 ))p;(s 1 - S2 )G ( x(gi) Xo r(S2) ) • 



We have shown earlier || that for (3 = 

A Q (X(tx),Y(t 2 )) = Ao(ti,t 2 ) =exp|^-jT t< 



G 



*oo-r(*) 



(is 



(C3) 



where t< = min(t 1 ,t 2 )- In computing the first order corrections in f3 for A ( g(X(t 1 ), K(£ 2 )), 
we shall proceed as in the previous appendix, by making an expansion 



ijm dp (s) 

~2 ds~ 



(C4) 



and taking the limit n — > oo, which allows us to make the replacement Po{s) — > 5(s) in all 
the integrals. The only term which requires a more careful treatment in the equation for 
the reduced influence functional is 



dAf,(X( Sl ),Y{s 2 )) ^ OAo(X( ai ),y( 8a )) = ^ 2 , 



9(si - s 2 ) 



9(si - s 2 ) 



(C5) 



which can thus be neglected. The reason is that Ao(X(si), Y(s 2 )) has a discontinuous partial 
derivative at s\ — s 2 = 0. Remembering that we need the influence functional for t\ = t 2 = t 
only and by using the obvious exact representation of the integral term 



dsi I ds 2 F(s 1 , s 2 ) 
o Jo 



t rsi ft rS2 

dsi ds 2 F(s 1 ,s 2 ) + / ds 2 ds 1 F(s ll s 2 ) 
o Jo Jo Jo 



(C6) 
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for an arbitrary function F(si,s 2 ) and by applying the rules described in the previous 
appendix we obtain that 



MX(t), Y(t)) = 1 + y f Q ds\ p {X{s), Y{s))G 

i[3T l rt 



'X(s)-Y( S y 



X 



+ 



4X Jo 



ds\ p (X(s),Y(s))[X(s) + Y(s)}G 



,( x( s )-Y( s y 

X n 



(C7) 



where G'(x) = dG(x)/dx. In the above evolution equation one can use either A (si,S2) or 
A.p(X(s) : Y(s)) on the r.h.s. An alternative way to derive this equation is to use the fact that 
for two given paths X(s 1 ) and Y(s 2 ) one has the explicit symmetry A (si,s 2 ) = A (s2? s i) 
and also that 



dP (s) dP (-s) 



ds 



ds 



(C8) 



before making the replacement Po(s) — > 5(s) as described earlier. 

This evolution equation can be easily solved as in the case of (3 = and the final answer 
for the influence functional is 



C(X(t),Y(t))=exp 
x exp 







/ ds 


H 







'X(s)-Y(s) 



- 1 



4X Jo 



ds[X(s) + Y(s)]G 



,fX(s)-Y(s) 



(C9) 



APPENDIX D: EFFECTIVE HAMILTONIAN AT FINITE TEMPERATURES 

Having obtained an expression for the influence functional at finite temperatures we 
can write down from the Feynman- Vernon path integral an expression for the effective 
Lagrangian: 



-iT l 



G 



X-Y 
Xo 



^ h lX + Y]G^ X - V 



4X 



X 



(Dl) 



We can introduce now the corresponding canonical conjugate momenta 



(D2) 
(D3) 



and construct the effective Hamiltonian according to usual rules 
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H = P X X + PyY 

1 



2M 
+ iT l 



Px _ fX ~ Y 



G 



4X 
X-Y 



X n 



U(X)- 



2M 



X 







-U(Y) 
(D4) 



The requantization of this effective Hamiltonian is straightforward (Px — * —ihdx and Py • 
—ihdy) and it is convenient to reorder the different terms in the Hamiltonian as follows: 



H 



(P x + P Y )(Px - Py) 



AMX Q 



a 



2M 
X 



+ U(X) - U(Y) 



Y 



X n 



(P x -P Y )+tT i 



G 



/X-Y 
V X 



- f 



(D5) 



Any other choice of ordering leads to an evolution equation for the density matrix p(X, Y, t) 
which does not conserve probability. One can consider alternative orderings, but in the final 
analysis these lead to slight renormalizations of various quantities, but to no qualitative 
effects. 
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Table I. The definition of our parameters as well as some of the limits used in this analysis. 



Parameters: 

Xq Characteristic scale over which Hi (X) changes. 

N Dimension of the Hilbert space of the intrinsic subsystem. We take the limit N ^> 1. 

k The bandwidth of the random matrix Hi(X). The average number of states coupled 
together at a given excitation energy is determined by the density of states: N ~ 
Kop(e). For k — > oo, we recover the full random matrix limit. 

Spreading width. For an initially uncorrelated state evolving under a random matrix, 
the average propagator decays as a(t) ~ exp(— TH/2h). 

(5 — 1/T Inverse thermodynamic temperature, defined through the density of states. 

p(e) The density of states. We use the form p(e) = p exp(/3e). When (3 = 0, the intrinsic 
system has a constant level density. 

G(x) Correlation function for the intrinsic states. It describes how far one must go in 
X before the intrinsic states are statistically uncorrelated. Typically one can use 
G(x) = exp(— x 2 /2), 1 — x 2 /2 or cosx. 

7 Friction coefficient obtained in the classical limit and in the full quantum dynamical 
solution, given by 7 = (3T l h/2MX^. 

D Diffusion constant, given by D = 2X 2 /(3 2 T l h. 
Limits: 



X — > 00 Weak coupling limit. 

Brownian motion limit. 



k/X -> 

Xo/rt -> 

7 = finite 



h — > Classical Limit. 

k q — > 00 Adiabatic Limit. This implies that the collective time scale X /V is much longer than 
the intrinsic one given by H/kq. 
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Figure Captions 



Figure 1. Time dependent flow S(t) associated with the time evolution of the density matrix 
in a linear potential, (top) The case of a Gaussian correlation G(x) = exp[— x 2 /2]. 
(bottom) The flow for a periodic correlation function G(x) = cos[x]. The stability of 
the lines of fixed points are indicated by the directions of the arrows. 

Figure 2. Fourier transform of the momentum distribution, d(s, 0, t), as a function of s, for F = 
and M = H = (3 = 1. From narrowest to widest we have X = 0.1, 0.5, 1, 2, oo. The 
Xq = oo (widest curve) is a Gaussian distribution. 

Figure 3. Fourier transform of the coordinate distribution d(0, k, t) for the case of a linear poten- 
tial. The solid curve is the result obtained from Eqs. |5.5| and |5.7| from integration to 



t = 100. This is compared to a Gaussian (dashes) obtained from the second cumulant 
alone. The top and bottom correspond to different correlation lengths Xq = 1 and 
X = l/10. 

Figure 4. Time dependent flow (S(t),K(t)) associated with the time evolution of the den- 
sity matrix in a harmonic oscillator potential, with the correlation function G(x) = 
exp[— x 2 /2]. The values of the parameters are to — 1 and 7 = 0.1 (top) and 7 = 0.5 
(bottom). 

Figure 5. Eigenvalue spectrum for the time-evolution operator of the density matrix p(r,s,t) 
for the harmonic oscillator with (3 > 0. For (3 = 0, the friction 7 vanishes, and all 
eigenvalues lie on the real axis. The A = solution corresponds to the equilibrium 
density matrix. 

Figure 6. Fourier transform of the momentum and coordinate distribution (they coincide in this 
case) for the harmonic oscillator for the case M = h = u = (3 = l here represented as a 
function of yl, where / is the action. The curves correspond to X = 0.5, 1, 2, 00 from 
the narrowest to the widest. The Gaussian distribution (Xq = 00), obtained from the 
Brownian motion limit in which only the quadratic cumulant is non-zero, is shown as 
a dashed curve. 

Figure 7. Time dependent flow associated with the time evolution of the density matrix in an in- 
verted harmonic oscillator potential, with the correlation function G(x) = exp(— x 2 /2). 
The parameters are to = 1/10 and 7 = 1/2. The flow of certain trajectories are shown 
with the arrows. 

Figure 8. Double well potential and eigenstates in the absence of the coupling to the intrinsic 
degrees of freedom. 

Figure 9. Strengths tunneled to the right hand well for = 0.5, 1 and 7r (upward from the 
lowest curve) and (3 = 0. 

Figure 10. Strengths tunneled to the right hand well for 1^ = it and (3 = 0,0.05,0.10,0.15,0.20 
(downward from the highest curve). 
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Figure 11. Density distributions at t=0 (the narrowest curve) and at time t=4 for Xo = 0.2, {3 = 
0, = 1: the continuous (dotted) curve corresponds to using the quadratic (Gaussian) 
correlation function. 

Figure 12. (Top) Temperature behavior of the complex eigenvalues of the double well potential 
with the smallest imaginary parts. (Bottom) Temperature dependence of the imag- 
inary part of the eigenvalues with smallest absolute value. One of the eigenvalues 
apparently vanishes in the limit T — > 0, which suggests that at T = dissipation 
decreases the tunnelling rate . One can see that the tunneling process is not always 
dominated by a single eigenvalue. 
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